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1 Introduction 


At the 17th J-PARC PAG, which was held on September 2013, we proposed 
the sterile neutrino search at J-PARC MLF [1]. After reviewing the proposal, PAG 
recommended to have a background measurement at the detector’s candidate site 
location in their report |2] to investigate whether the background rates can be man¬ 
ageable for the real experiment or not. Therefore, we have performed the background 
measurements (MLF; 2013BU1301 test experiment) during the summer of 2014, also 
following the 18th J-PARC PAG recommendations [3] , and the measurements results 
are described here. 


2 Background Measurement at the MLF 3rd Floor 


2.1 Goal and Strategy 

The goal of the 2013BU1013 measurements is to check the background rate, es¬ 
pecially related to the beam, written in the P56 proposal [1] . Table is the reprint 
background summary table from that proposal. Here we assumed that the two detec¬ 
tors, with 25 tons liquid scintillator hducial volume (50 tons in total), are set on the 
MLF 3rd floor with the baseline of 17 meters. We also applied the selection criteria 
shown in Table (Reprint of Table 7 in m) to select the Inverse Beta Decay (IBD) 
signal events efficiently, which was caused by the oscillated signal from the muon 
decay at rest. { ^ + Ue + i7fj. (at the MLF mercury target); i7^ —>■ (oscillation 

during 17 m)) —>■ he; ry + P —t e+ J- n (at the liquid scintillator detector). The back¬ 
ground rates from “beam associated fast neutron” and “Accidental events” are to be 
checked. These numbers were estimated by the MLF 1st floor measurements and the 
MG simulation in the proposal, therefore the direct measurement was recommended 
by the PAG. 

The strategy to check the numbers are; 

1. Background rates are measured by a 500 kg plastic scintillator and other small 
scale (less than 10 kg) size detectors described in the appendix. 

2. The rates are scaled to the P56 detector with the hducial mass of 25 tons at 
hrst, then multiplied twice for the two detectors. The IBD selection criteria are 
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Table 1: A reprint table from the P56 proposal jTj. Numbers of events of the 
signal and backgrounds with total hducial mass of 50 tons. 


Cut Condition 

Cut Efficiency 

1.0^ ^lO^S 

74% 

^A^delayed <12MeV 

78% 

20<Eprompt ASOMeV 

92% 

^^delayed ^100/XS 

93% 

AVTX<60cm 

96% 

Total 

48% 


Table 2: A reprint; IBD Selection criteria and efficiencies for the oscillated signals 
for high Am^ case. 


also applied. The scale is based on the MC simulation since the 25 tons detector 
has self-shielding effects against gammas and neutrons, and the volume ratio or 
the surface ratio between the 500 kg (or small size) detectors and the 25 tons 
detector does not account for the rate alone. 

3. In this report, we assume that the 25 ton detector has shielding effects only 
from the stainless tank and the liquid scintillator inside the tank. There are no 
additional shields made from iron, concrete or mineral oils. 

Figure [T] shows the current 25 tons detector design described in the Appendix]^ 
Before backgrounds’ arriving at the acrylic vessel, there is a buffer region made by 
the liquid scintillator, whose thickness is 50 cm. This thickness is almost the same as 
the height and width of the 500 kg detector. 

^ As describe in the proposal, the 3rd floor of the MLF is the maintenance area for such as a 
mercury target, beam monitor equipments, and so forth. Therefore, the detector should be moved 
to outide of the building for the maintenace works, which is held at least once per year. The detector 
size is restricted by the entrance of the building. 
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Figure 1: Current design of the P56 25 tons detector. Left; schematic view, right; 
3D drawings. The “n -7 shield” in the left cartoon is not considered in this report, 
thus only stainless tank and scintillator inside the tank are assumed to shield the 
neutrons and gammas. No additional shields such as irons, concrete or mineral oils 
are assumed. 

2.2 Measurement Points and Operation Time 

Figure 1^ shows the overview of the MLF building (left), and the measurement 
points of the 2013BU1301 experiment (right). The 500 kg plastic scintillator detector 
and their calibrations are described in Appendix A. The baselines from the mercury 
target are 17m, 20m, 34m for Pointl, 2 and 3 in Fig. respectively. We performed 
the background rate measurement for two weeks per each point. 



Neutron Hgtar^^ 
(& neu^ino source) 


Pointl 


Points 


3GeV p beam 


Point2 


Figure 2: A schematic view of the the MLF facility in J-PARC (left), and the 
measurement points (red circles) of the 2013BU1301 experiment, (right; those are 
written as “Pointl”, “Point2” and “Point3”)) 

Only the result of “Point2” is described here since the amount of the background 
rate at “Pointl” is much larger than the other points due to the proton beam multiple 
scattering at the muon target of MLF [1] , and the baseline for “Point3” takes longer 
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time for the sterile neutrino search in the current conhguration using the 25 tons 
detectors due to the statistics. 

In order to extrapolate the background rates from the test experiment to those 
of the proposed P56 25 tons detector, we assume 24 meters baseline for the 25 tons 
detector in this report since the realistic two detectors are to be put around Point2 
to manage the background rates. However, we cannot put the detectors in the area 
of the beam upstream of the Point2 since there are many equipments for the MLF 
facility there. The possible detector location is being discussed in details 
with MLF facility people, and will be determined with all constraints. 

During the test experiment, we found that the default operation time of the MLF 
facility is 5000 hours / year. Therefore we changed to the operation time from 4000 
hours written in the proposal to 5000 hours / year for the sensitivity calculation. 


2.3 Background from Beam Fast Neutrons 

One of the main purposes of this background measurement is to directly measure 
the Michel electron background induced by beam fast neutrons, which was indicated 
by the previous background measurement at BL13 [T]. According to the Geant4 
Monte Carlo simulation, fast neutrons whose kinetic energy are larger than 200 MeV 
can produce charged pions. The flux of such fast neutrons at MLF 3rd floor is 
estimated to be some order of magnitude smaller than that at BL13. We measured and 
conhrmed that the beam Michel electron background is low enough for the experiment. 
We will briefly describe the basic idea and backgrounds of this measurement. 

Figure 1^ shows the dehnitions of the “signal” and “background” of this measure¬ 
ment. The signal is a Michel electron induced by beam fast neutrons. Fast neutrons 
coming on the (proton) bunch timing hit our detector and produce pions. These pions 
then decay into Michel electrons (n + p{orC) X + tt’*', then tt/ r+ —>■ e"*"). The 
signature of the signal is thus the coincidence between a neutral activity on the bunch 
timing and a “prompt signal” about 1 /is later from the beam timing. Backgrounds 
of this measurement are clipping cosmic muons, Michel electrons from cosmic muons 
and neutral particles (gammas and neutrons) from cosmic rays as was shown in Fig.[^ 
and they come either during beam-on or beam-off. On the other hand, signal comes 
only when the beam is on. The basic idea of this measurement is thus to extract 
signals from backgrounds by subtracting beam-off activities from beam-on activities. 

Based on the concept of the measurement, we took the following data set: 

• beam-on: To observe activities on the bunch and just after bunch timing (the 
timing dehnition is described later); 

• 20ms-later : To subtract backgrounds from beam-on data 

We took data with a 20 ms delay after each beam bunch spill. Because of the 
PMT gains, efficiency of veto counters and other detector response are exactly 
the same with the last beam spill, we can extract the backgrounds from the 
beam signals with less systematic uncertainties; 

• beam-off: To evaluate truly beam unrelated background; 
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Figure 3; The definition of “signal” and “backgrounds” of this measurement. Signal 
is a Michel electron induced by beam fast neutrons. Backgrounds are clipping cosmic 
muons, Michel electrons from cosmic muons and neutral particles from cosmic rays. 
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• cosmic muons: To calibrate the detector. 


2.3.1 Measurement 

The search for beam neutron induced Michel electrons was performed by detecting 
their prompt signals. Figure shows the correlation between the energy and the 
timing of the events observed at Point 2. In this measurement the prompt signal is 
defined by: 

• 20 < i?[MeV] < 60 

• 1.75 < t[ns] < 4.65 from the rising edge of the first beam bunch. 

As described before, we compare the beam-on and also 20ms-later data to subtract 
other activities with less systematic uncertainties. A huge number of the clipping 
muon background is rejected by applying the charged veto cut. Figure shows the 
energy distributions of events in the prompt timing window, 1.75 < t[/is] < 4.65 from 
the beam bunches, and 20 ms later, before and after applying the charged veto cut. 
Without applying the charged veto cut, we observed (1.68±0.03) x 10“'^/spill for beam 
data, and (1.64T0.03) x 10“^/spill for 20ms-later data. Applying the charged veto cut, 
we observed (1.58 ± 0.09) x 10“®/spill for beam data, and (1.52 ± 0.09) x 10“®/spill 
for 20ms-later data. The numbers of events in the prompt energy range are both 
consistent between beam and 20ms-later data either with or without applying the 
charged veto cut. 

To improve the sensitivity, an additional cut was applied before the subtraction. 
Figure 1^ shows the energy distribution on the bunch timing when fast neutrons pro¬ 
duce charged pions in the 500 kg detector (MC). We can observe some events on 
the bunch timing associated with beam Michel electron backgrounds. On the other 
hand, as was shown in Fig. most of the beam unrelated backgrounds do not have 
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Figure 4: Correlation between energy and timing of the events observed at Point 2. 
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Figure 5: Energy distributions of events in the prompt timing window and 20 ms 
later, before and after applying the charged veto cut. 
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Figure 6: Energy deposit on the bunch timing for beam Michel electron with the 
neutron kinetic energy 300 to 500 MeV (flat). Note this is MC. Red, blue and black 
lines are the energy deposit from neutron interactions, pion decays and their sum in 
each event. A peak around 3 MeV in blue histogram corresponds to the muon kinetic 
energy, 4.2 MeV, from stopped pion decay including the Birks’ quenching. Birks’ 
quenching, light attenuation in the scintillator and other detector responses such as 
resolutions and threshold effects wer e imp lemented to the Geant4 based Monte Carlo 
simulation as described in Appendix |A.4 


any activities on the bunch timing. Michel electrons from cosmic muons can have 
an activity on the bunch timing only when the parent muon comes on the bunch 
timing accidentally. However, since the muon is a charged particle, it can be easily 
rejected by external veto counters. We can thus strongly suppress backgrounds by 
requiring on-bunch activities without hits in veto counters. At least one on-bunch 
activity (Edep>4 MeV) without hits in veto counters was required. Figure shows 
the estimated selection efficiency of this on-bunch cut as a function of the incident 
neutron kinetic energy based on MC. Though most of the events have more than 4 
MeV energy deposit at on the bunch timing, a part of the events are rejected by 
self-vetoing. The selection efficiency has slight dependence on the incident neutron 
kinetic energy. Because we do not know the energy spectrum of the incident neutron 
well, we assumed the selection efficiency, Conbunch = 0.9. 

Figure [^shows the energy distributions after applying the cut (the 500 kg detector 
data). The beam unrelated distribution was obtained with the energy distribution of 
20ms-later data without applying the on-bunch cut, and the accidental coincidence 
probability: hit rate of neutral activities on the bunch timing. The mean hit rate 
without veto activities on the bunch timing was 3.1%. The observed event rate during 
beam-on was (4.60 ± 1.53) x 10“^/spill, while the estimated event rate by beam 
unrelated activities was (4.91 ± 0.28) x 10“^/spill, and both rates are consistent. By 
considering the efficiency of the on-bunch cut, Conbunch = 0.9, the upper limit of the 
event rate of the beam Michel electron is thus 2.5 x 10 ^/spill (90% C.L.). Because 
we set the timing window, 1.75 < t[/is] < 4.65 from the beam, for the beam Michel 
electron, the obtained upper limit is thus equivalent to 4.7 x 10“^/spill (90% C.L.) 
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Figure 7: Estimated selection efficiency of the on-bunch cut as a function of the 
incident neutron kinetic energy. 


for the timing range, 1 < t[/is] < 10 from the beam. 

We then extrapolate the obtained upper limit to the 25 tons detector. By consider¬ 
ing the difference of the beam power (300kW —IMW) and the detector acceptanc^ 
the obtained upper limit is equivalent to 1.1 x 10“®/spill/detector/MW (90% C.L.). 
The upper limit is about 53 times larger than the event rate assumed in the proposal, 
and it is simply limited by the target mass and the measurement period. 


2.4 Background from Accidental Coincidence 

Another important background is due to accidental coincidences. In order to es¬ 
timate the accidental background, an absolute rate of the background measurements 
on the prompt and the delayed region with the 500 kg detector are crucial. Measure¬ 
ments with small size detectors are also important to estimate the contents (BID) of 
the prompt background and to estimate the rejection ability of the lead blocks of low 
energy gammas. 

We first explain the background measurements for the prompt signal region, and 
then for the the delayed region. 


2.4.1 Background Measurements for Prompt Backgronnd Region 

The estimation of the contents and rate of the prompt region from a small size 

^ We assumed fast neutrons, which produce beam Michel electrons, are directly coming from the 
mercury target. The lead shield below the 25 tons detector, described in Section |2.4.2[ was also 
considered. The scale factor by the detector acceptance is estimated to be 7.6 in average over two 
detectors. 
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Figure 8: Energy distributions after applying the on-bunch cut. The beam unrelated 
distribution was obtained with the accidental coincidence probability, 3.1%, and the 
energy distribution of 20ms-later data without applying the on-bunch cut. 


liquid scintillator and Nal(Tl) detectors is shown first, then the measurement with 
the 500 kg detector for the prompt region is also described. The results of the small 
size detectors and 500 kg detector are scaled to the proposed 25 tons detector at the 
end. 


PID and Energy Measurements with Small Size Detectors 


Appendix [^ shows the complete description of the measurements with the small 
size detectors at Tohoku University, therefore only a summary of the measurements 
is shown here. 

A Sodium Iodide (Nal) scintillator counter (a cylinder with 2” of diameter and 
2” height) was hrst used to measure the gamma flux that contributes to the prompt 
event of the accidental background. Cosmic ray muons are rejected using plastic scin¬ 
tillators, which surrounds Nal counters, with an efficiency better than 99%. Figure]^ 
shows the results of the measurement, where the remaining gammas energy spectrum 
can be described by two exponential functions with decay constants of 3 and 26 MeV 
(the gammas for MC were generated following the two exponential functions from 
the Nal surface assuming an isotropic direction). Comparing data and MC spectra, 
these components rates were measured to be 150 and 25 Hz/m^, respectively. 

Another important measurements has been done using a small liquid scintillator 
(NE213) which can separate neutrons from gammas efficiently. A cylindrical alu¬ 
minium housing (5” of diameter and 2” height), with white painted inner walls, was 
hlled with NE213, closed with a glass plate and attached to a 5” PMT (R1250-03). 


The NE213 detector is also surrounded by the plastic scintillators, and Fig. 10 shows 
the results of the measurement. 

With this set-up, the fast neutron flux detected above 20 MeV is of (1.28±0.05) x 
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Figure 9: Nal energy spectrum for the remaining events after applying the veto cut. 



Energy (MeV) 


Figure 10: Reconstructed energy distribution for events with low or none energy 
deposition in the veto. The data and MC components of neutrons and gammas, 
selected as explained in Section B, are also compared. Data events without veto 
applied are also shown for comparison. 
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10“^ Hz (statistical uncertainty), while the MC gives a rate of 1.12 x 10“^ Hz. These 
numbers are close to the one used in the proposal pQ. For gammas, in the same energy 
range, the measured rate is of (1.18 ±0.04) x 10“^ Hz (statistical uncertainty), while 
the MC gives 0.95 x 10“^ Hz. Therefore, the data and MC above 20 MeV for both 
gammas and neutrons agrees within 20% of the uncertainty. For the MC presented 
here and for the one in the hrst proposal, the same cosmic neutrons generator was 
used. The generator’s flux and spectrum were tuned with the Tohoku University’s 
Reactor Monitor detector, composed of 200 litters of liquid scintillator [I] . 


Measurement for Prompt BKG with 500 kg Detector 


Figure shows the energy distribution of beam-off data after applying charged 
veto and the Michel electron cuts. And the Monte Carlo estimation with cuts is 
overlaid. The cosmic muons are eliminated by applying the charged veto cut because 
of the high veto efficiency (>99.5%) as described in Appendix A.5 The Michel elec¬ 


trons can be rejected by detecting the parent muons coming earlier than the prompt 
events, and the remaining events are composed of cosmic gammas and neutrons. 

The 25 tons detector for the sterile neutrino search will have PID capability, which 
can reject fast neutrons relative to electrons and gammas by a factor of more than 
100 pp. Neutrons in the remained events are thus not harmful and only gammas can be 
a prompt background. The flux ratio and energy distributions of the cosmic gamma 
and neutron were measured by using small size liquid scintillator and the Nal(Tl) 
detector at Tohoku university as described in the previous section. By using the y/n 
ratio, the number of gammas in the remained events is equivalent to (6.4±0.5) x 10“® 
/2.5/is. 



energy[MeV] 


Figure 11: Energy distribution of beam-off data after applying charged veto cut and 
the cosmic Michel electron cut. 


Note that the predicted rate from the measurement at the RCNS (Tohoku Univer- 
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sity) is consistent with the rate measured by the 500 kg detector at the MLF within 
6%. Thus, a reliable extrapolation to the proposed 25 tons detector is possible, since 
the MC model agrees for different detector types (Nal, NE213, Reactor Monitor, and 
500 Kg plastic scintillator) at different locations (RCNS and MLF). 


Background at the proposed P56 25 tons detector 

As a hnal step, the background measurements explained in the previous sections 
were scale to this 25 tons proposed detector. When generating high energy gammas, 
following the measured spectrum and the same MC generators for the 500 kg detector 
for the proposed detector, the rate of events remaining after the selection criteria is 
of 3.8 X 10 events/spill/9/is, a value 29 times higher than the last proposal. This is 
simply because we didn’t consider this contribution. Figure 1^ shows the spectrum 
of such events. The estimated number of neutrons events, using the same rate as in 



Figure 12: High energy MC gammas energy distribution for the proposed detector. 

the previous proposal, is 32.2 in 4 years, considering 4000 hours of operation in each 
year. 


2.4.2 Background measurements for Delayed Signal 


In this section we report on background measurements and background estimation 


of the 25 tons detector for delayed signal. As seen in section 2.4.1 the background 
rate for the prompt signal is about 30 times higher than the last proposal. Therefore 
we have to manage the accidental background rate reducing the background rate for 
the delayed signal to a value more than 30 times smaller than the proposal. The 
schematics of the 25 tons detector is shown in Fig. Outside the acrylic vessel is 
the buffer region, which is separated into two layers with a black sheet. The most 
outside region is used as a veto cut. 

The energy region for the delayed signal is 7 < i7[MeV] < 12. This range is 
different from the region taken in the previous proposal to reduce the amount of 
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background. The effect of tightening energy selection is discussed in Section 3.1.1. 

Backgrounds for the delayed signal were also measured with the 500 kg scintillator 
detector at the 3rd floor in MLF. There are two types of backgrounds for the delayed 
signal. 

• Gammas by neutron captures in the materials outside of the detector. The 
gammas directly hit the detector from outside and mimic a delayed signal. 

• Neutrons which come into the detector from outside, are captured in the detec¬ 
tor and emit gammas of which energy is same as delayed signal. 

Based on the measurement, we estimated MC models of backgrounds for gammas 
and neutrons. Using the MC models, we estimated amounts of background in the 25 
tons detector. 


Delayed Background from Beam Gammas 


As a preliminary measurement we observed an unexpected background rate around 
Point 2 and 3 with a small scintillator and an oscilloscope. Around Point 2, the back¬ 
ground rate was more than 1 kHz with a threshold equivalent to ~ 1 MeV. This rate 
was more than 10 times higher than that of other locations including Point 3. The 
location and result are shown in the top part of Fig. Because this rate decreased 


to 200 Hz when 5 cm-thick lead blocks were put below the small scintillator, this 
background can be regarded as gammas. In addition, this rate decreased to 100 Hz 
with 10 cm-thick lead blocks under the small scintillator. 

On the other hand, the rate decreased to 150 Hz when 5-cm thick lead blocks 
surrounded almost all over the small scintillator. Because lead blocks of bottom side 
had larger shield effect than those of other side, gamma background was coming from 
the floor. (See also Appendix C for more details.) 

Measurements with the 500 kg detector also indicated the background rate at 
Point 2 was several ten times higher than the rate at Point 3 in the energy region 
of 7 < i7[MeV] < 12. Therefore it is important to understand the source of gamma 
background at Point 2 and how to suppress it. 

The energy spectrum for the delayed signal at Point 2 is shown in the black dots 


of Fig. 14, In the energy region of 7 < i7[MeV] < 12, the background rate with the 
500 kg detector is (2.06 ± 0.01) x 10“^/2.5/is/0.3MW. To explain this background, 
we estimated a MC model. 

As shown in Fig. the energy distribution does not change till 40 ms after the 
beam bunch. Therefore the background rate for the 2.5 ps time window is linearly 
scaled to the rate for the 100 ps time window. 


MC model for gamma background 

In order to construct a reliable MC model for gammas production, the followings 
are assumed. As shown in the bottom part of Fig. 13, the remote maintenance space 
of the mercury target exists under Point 2 and 3. Neutrons generated by the beam are 
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Figure 13: Top part of the figure shows the top view of measurement location and the 
result of preliminary measurement. In the cyan area around Point 2, the background 
rate was more than 1 kHz. In the gray area, the background rate was about 100 Hz. 
Bottom part of the hgure is the side view of measurement location. The floor of Point 
2 is made of concrete which has 120 cm thickness and is thinner than that of Point 3 
by 30 cm. 



energy[MeV] 


Figure 14: Energy distribution of background for delayed signal and other activities 
measured with the 500 kg detector 
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Figure 15: Energy distribution of gammas in the MC model of Point 2 


captured at the ceiling of the room which is made of concrete. Gammas are emitted 
as a secondary product of the neutron capture by the concrete, and consequently 
they reach Point 2 and 3. Because the concrete thickness at Point 2 is about 120 cm, 
which is removable for maintenance and thinner than at Point 3 by 30 cm. Point 2 
have a higher background rate. 

Having these assumptions in mind, the energy and angle distributions of gamma 
of the MC model are shown in Fig. 15 and 16 Gammas reaching Point 2 have energies 


reflecting the neutron capture by concrete and vertical gammas reach the detector 
easier than the ones with large angle. 

In the MC simulation, these gammas are generated uniformly on a floor around 
the detector. We validated this MC model, comparing between measurements and MC 
with the small plastic scintillator in various shield configurations (see Appendix |C.4[ ). 
We also validated the MC model with data taken by the 500 kg detector. As shown 
in Fig. 14, the energy spectrum is well explained by the MC model (blue area). 


environmental gammas (green area) and cosmic rays (red area) 


Extrapolation to the 25 tons detector for gamma backgronnd 


Using the MC model, we estimated the background rate of the 25 tons detector. 
The background rate, A^expected) is 


expected 


^MC ArSOOkg 
■‘'obs ’ 


/ySOOkg 


( 1 ) 


where A^mc i® calculated number of events of the 25 tons detector with MC, 
is the calculated number of events of the 500kg detector with MC and is the 

observed background rate. Without any shield, A"gxpg(,^gj is 1.5 /100/is/MW/detector 
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cos 6 


Figure 16: Angle distribution of gammas in the MC model of Point 2. cos 6 = 1 
means that gamma with a vertical direction. 


in the energy region of 7 < i7[MeV] < 12. This value is 32 times larger than the 
value of the proposal. To compensate this increment, at least ~ 6 cm thick lead 
blocks are needed. As conhrmed in Appendix C.4[ adding a lead shield below the 
detector is effective to reduce this gamma background. Note that we tested the effects 
of 10 cm thick lead blocks in the conhgurations. In this report, we assume to put 
12.5 cm thick lead shield under the 25 tons detector to suppress the background, and 
to compensate the prompt background increase as well. Thus, N^^^^cted reduced to 


1.2 xlO ^ /100/rs/MW/detector. 1^ 


Delayed Background from Beam Neutrons 

As shown in Fig. the 25 tons detector have a 50 cm thick buffer region that 
acts as a self-shield. This shield effect for the low energy neutrons (less than 10 
MeV) is shown in Fig. (See Sec. 5.8 of [T]) and its neutron rejection power is 
shown in Table including the energies higher than 10 MeV. Low energy neutrons 
(kinetic energy En < 10 MeV) is signihcantly reduced by self-shield of the 25 tons 
detector. On the other hand, high energy {En > 10 MeV) neutrons can reach the 
hducial volume and are thermalized and captured by the Gd. Therefore neutrons 
with En > 10 MeV become background for delayed signal. 

Because such neutrons come at beam on-bunch timing only, we estimated the 
background rate of neutron at the 25 tons detector from on-bunch energy spectrum 
measured with the 500 kg detector. We use the energy spectrum above 15 MeV since 
the gammas are dominant below 15 MeV. For on-bunch timing, we use 3 /is time 

^ As discussed later, we apply AyTXoB-deiayed cut to reduce neutron background. The back¬ 
ground rate of beam gammas is also affected by the cut and reduced by 2.0% as well as neutrino 
signal. However, the effect of the reduction can not be seen in this digits. 
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Figure 17: 


Self shield effect of liquid scintillator as function of the shield thickness. 


Initial energy [MeV] 

0-20 

20-40 

40-60 

60-80 

80-100 

rejection power 

75 

13 

6.1 

4.4 

3.7 

Initial energy [MeV] 

100-120 

120-140 

140-160 

160-180 

180-200 

rejection power 

3.3 

3.1 

2.9 

2.7 

2.6 


Table 3: Rejection powers of selfshield of the 25 tons detector. Rejection power is 
ratio of the number of incoming neutrons at the surface of the 25 tons detector to the 
number of neutrons reaching the surface of the hducial volume. 


window after the rising edge of the hrst beam bunch. In Fig. 2D plot for energy 
and timing of activities measured with the 500 kg detector is shown. Energy of on- 
bunch activities is projected to the black cross in Fig. 1^ An estimating method is 
as follows; 


1. We estimated a neutron flux which explains the measurement of the 500 kg 
detector with MC. 

2. Amount of activities in the 25 tons detector was estimated using the flux and 
MC. 


For the following estimation, we assumed that all activities of on-bunch timing in 
the 500 kg detector are caused by neutrons. This assumption could be conservative 
to estimate amount of neutron background because there could be some fraction of 
gammas as well. 
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Figure 18: Energy distribution of on-bunch timing with the 500 kg detector. The 
black cross is the measured energy spectrum and the red one is a MC fitted to the 
measured one. 


Neutron flux 


To explain the measured spectrum, we estimated an energy spectrum of the neu¬ 
tron flux. The geometry used is shown in Fig. [T^ We assumed that neutrons come 
directly from the mercury target. The flux, (j){E^), is expressed as 


a 


(j){En) = — exp(-F;n/30MeV), 


( 2 ) 


where En is kinetic energy of neutron and a is a number of neutrons that are generated 
at the mercury target in a spill, a was estimated to be 387/spill as a ht result. To 
calculate the visible energy spectrum of the MC, we considered the quenching effect 
of plastic scintillator using the Birks’ law as well as the neutron reactions inside 
the scintillators in GEANT4. As shown in Fig. 18, this flux reproduces well the 
measurement above 15 MeV. 


Estimation of neutron background for the 25 tons detector 


Rate calculation 

Using the estimated flux, we calculated the neutron rate, which creates the fake 
delayed activity of the 25 tons detectoi]^ This fake delayed gammas from neutron- 
captures mimic the neutrino signals if prompt fake gammas also exist accidentally. 

® Note that the Birks’ law as well as the neutron reactions inside the liquid scintillator in GEANT4 
are also included to calculate this rate. 
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Figure 19: Geometry of the mercury target and the 500 kg detector 


The geometry used is shown in Fig. In the MC, the 25 tons detector is located 
above a 12.5 cm thick lead shield. Such a conhguration is used for the MC gamma 
background as well. 

As a result of the calculation, we have a rate of the background in delayed region 
of 0.0163/spill/MW/detector. 


Further Reduction Methods 


This rate is much larger than the background from beam gammas, thus we have 
to consider further reduction methods. 

Fortunately, 0.0155/spill/MW/detector out of 0.0163 of the background has hits 
which have more than 0.5 MeV energy inside the black-sheet region in Fig. on 
the beam bunch timing due to the proton recoils of the beam neutrons. This re¬ 
coil proton hit position on the bunch timing and the delayed fake gamma position 
after the neutron thermalization has a strong spatial correlation. (The red line in 


Fig. 21 is the distance between the hit position of the recoil proton and that of the 
delayed gamma). This correlation (we call this variable Al/TXoB-deiayed) reduces the 
background from neutron-capture by 97.4% when setting 110 cm as the cut value. 
Therefore the remained background out of 0.0155/spill/MW/detector is 0.0155x0.026 
= 0.0004/spill/MW/detector. 

There are invisible energy {E < 0.5 MeV) hit events on the beam bunch timing 
inside the black-sheet region (0.00077/spill/MW/detector). However, most (97%) of 
them have hits in the veto (outside the black-sheet) region {Eveto > 0.5 MeV). 
Therefore, we can cut them using the information. The remained background rate 
out of 0.00077/spill/MW/detector is negligible compared to the former one. 

Using these two further reductions, the hnal delayed background rate from the 
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Detector 



Figure 20: Geometry of the mercury target and the 25 tons detectors 


beam neutrons is 0.0004/spill/MW/detector. 

3 Realistic Rate Estimation for the P56 Experi¬ 
ment using Data 


3.1 Additional Selection Cuts from the P56 Proposal 


3.1.1 Tightening Energy Selection for Delayed Activity 

By the measurements at MLF 3rd floor, beam-induced gammas from the floor 
were observed, and they contribute much to accidental events. To reduce the beam- 
induced gammas and to improve the S/N ratio, 6 MeV of lower energy cut for delayed 
events in the proposal is changed to 7 MeV in this analysis. Figure. 22 shows delayed 


energy spectra of the beam-induced gammas (red), gammas from thermal neutron 
capture on Gd (black), and e+ from €.~Y‘^Ng,s. reaction (blue) in the left hgure, 

and their cut efficiency curves in the right plot. Applying 7 MeV as the lower energy 
cut, 73% of the beam-induced gammas can be rejected while the signal efficiency can 
be kept to be 91 % of the 6 MeV cut. The signal cut efficiency for delayed events is 
changed from 78% to 71%. In the case of the ^^G(r'e, e“)^^Ag.s, reaction, the delayed 
energy cut efficiency is reduced to 85% of the 6 MeV lower cut. 


3.1.2 Additional Lifetime Cuts 

According to the 2013BU1301 test experiment, the background for the prompt 
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Figure 21: AVTX distributions of on-bunch and delayed signals (MC); the red line 
shows the distance between on-bunch hit and delayed hit positions from on-bunch 
neutrons while the blue line shows the distance between on-bunch and delayed hits 
from the uncorrelated events. 


signal region is dominated by gammas induced by cosmic rays. Therefore, the time 
structure of the prompt background is flat (Fig. |^. On the other hand, those for 
the delayed signal region is dominated by the beam induced gammas and neutrons 
as described before. For gammas, the measured time structure is also flat, while that 
for the beam induced neutrons has the same lifetime as neutrino induced events since 
beam neutrons come to the detector at the proton beam timing and are predominantly 
thermalized. 

At the P56 proposal [1], we applied the lifetime cuts for the prompt ( 1 < Atp < 
10/is) and the delayed activities {Atp < Atd < 100/is) separately. We call this as 
“baseline cuts” hereafter in this subsection. The lifetime of the oscillated signal hits 
are shown in Fig. [2^a) (prompt) and (b)(delayed). However, the tightened cut can 
be added to the baseline cuts using two dimensional correlation between the prompt 
and the delayed hits time from the proton beam. Table 0 and Fig. gc),(d),(e) 
show the features of the lifetime of events in three categories, the neutrino signal, 
the accidental from gammas (prompt) -|- gammas (delayed) and those from gammas 
(prompt) -|- neutrons (delayed). 



Prompt 

Delayed 

Comments 

Neutrino Events 
Accidental! 
Accidental 2 

Exponential (2.2/is) 
Flat 

Flat 

Exponential (30/rs) 
Flat 

Exponential (30/rs) 

7 -I -7 case 
7 -I- n case 


Table 4: Features of the hit time distribution for three categories. 


In order to maximize the rejection factor of the accidental backgrounds and to 
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Figure 22: Delayed energy spectra of beam-induced gammas and gammas from 
thermal neutron capture on Gd(left) and their cut efficiency curves (right). Blue 
dotted lines show the new lower cut. 


minimize the loss of the signal efficiency, the likelihood function based on the lifetime 
of the signal hits is created as follows; 


This likelihood is calculated in each event using hits time of the prompt and the 
delayed activities. Probability Density Functions (PDFs) to calculate the probability 
for the prompt and delayed activities are shown in Fig. [2^a) and (b). (note that 
these distributions are still not applied for baseline cuts). 

The calculated likelihood distributions, and the signal efficiency of these events 
(three categories) are shown in Fig. 24 As shown, we can reject more than half of 


the 7 -l- 7 type accidental background, while 91% of the signal is kept, applying a 
likelihood cut of 11.0. 

We set the cut value of the likelihood at 11.0 in this report. 


3.1.3 ADTWnR— 


OB —delayed 


Cut 


In section 2.4.2, we already discussed the effects on the reduction of the neutron 


backgrounds using the ADTAoB-deiayed cut. In this subsection, the topic is totally 
different from that section. The topic is not the background rate but the signal 
inefficiency for the oscillated neutrino signals due to the ADTAoB-deiayed cut. 

Although the ADTAoB-deiayed cut is a quite effective to reduce the neutron back¬ 
ground, some of the oscillated neutrino events are also cut accidentally. In short, if 
the space in a spill are overlapped between oscillated neutrino signal and background 
neutron events accidentally, the signal is cut. We have to consider two following cases; 


■‘on-bunch -|- delayed neutron events’' 


As mentioned in the section 2.4.2,0.0155/spill/MW/detector has both on-bunch 
and delayed activities from beam neutrons. This means that 1.55% of the 
oscillated neutrino signal events have the additional on-bunch -I- delayed neutron 
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Figure 23: (a) Hit time distribution (toy MC) of the oscillated signal (prompt), 

(b) Hit time distribution (toy MC) of the oscillated signal (delayed), (c) Correlation 
between prompt and delayed hits for the signal events after the baseline cuts, (d) Cor¬ 
relation between prompt and delayed hits for the accidental background (7+7) events, 
(d) Correlation between prompt and delayed hits for the accidental background (7 -|- 
neutron) events. 


activities as shown in the bottom of Fig. 25 If the vertex position of the 


neutron events are overlapped with the signal neutrino events, AHTXoB-deiayed 
cut rejects the signal. The blue line in Fig. 21 shows the AHTAoB-deiayed 
distribution for uncorrelated events such as this case, and an accidental cut 
efficiency (AHTAon-deiayed < HO cm) for the uncorrelated events are 11%. 
Thus, the signal inefficiency in the case is 0.0155 x 0.11 = 0.17%. 

2. “on-bunch only nentron events”; 

In the case that beam neutrons create the hits on the beam bunch timing but 
do not make any delayed gammas (e.g. due to escaping the detector region), it 


also produces the signal inefficiency. (The top schematic in Fig. 25). According 


to the rate calculation, 16.3% of the signal neutrino events have such additional 
“on-bunch only neutron”, therefore we have a signal inefficiency of 0.163 x 0.11 
= 1 . 


In total, we have 2.0% for the signal inefficiency. 
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Figure 24: The calculated likelihood distributions (left) and the efficiencies of events 
in three categories as a function of the likelihood cut value (right). 

3.2 Grand Summary of the Background at P56 

3.2.1 Summary of Selection Criteria 

Table shows a summary of the new selection criteria and the cut efficiency for 
the oscillated signals with IBD (high Am^). The lower energy cut for delayed events 
is changed, and two new cuts (the AVTXos-de^at/ed cut and life time cut) are added 
to the selection criteria in the proposal. This is because the current estimation of 
the accidental prompt events have a 30 times larger background rate than that in 
the proposal, based on the observation in this background measurement at MLF 3rd 
floor. Then, the total selection efficiency for the oscillated signal in the high Am^ 
region is reduced to 38% from 48% in the proposal (19% reduction from the proposal). 
However, the accidental background can be much reduced to an amount comparable 
with the proposal. 


Cut Condition 

Proposal 

New 

^^prompt ^10/iS 

74% 

same 

Q^^delayed <12MeV 

78% 

71%(Ed>7MeV) 

20<Eprompt <60MeV 

92% 

same 

^^delayed ^100/iS 

93% 

same 

^^T^^prompt—delayed ^60cni 

96% 

same 

AVTXqb- delayed >110cm 

- 

98% 

Life Time<ll 

- 

91% 

Total 

48% 

38% 


Table 5: Summary of selection criteria and efficiency of the oscillated signals for 
high Am^ region. 
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On-bunch neutron event cases 


On-bunch only neutron event 

On bunch Prompt Delayed 


1 


Neutron in 
Black sheet 


I r 

■ Inverse beta ■ 

I decay signal I 


Time 


Signal inefficiency 
due to AVTXoB-delayed CUt 


0.163x0.11=1.8% 


On-bunch + delayed neutron event 

On bunch Prompt Delayed 

“f T rr ^ 

■ Neutron in ■ Inverse beta ■ fon bunch 
f Black sheet I decay signal I i related 


0.0155x0.11=0.17% 


Total 


1.97% 


Figure 25: The cases signal is rejected due to Al/TXoB-deiayed cut and its inefficiency. 










3.2.2 Accidental background Summary 

Table shows an accidental background rates summary per detector. Concerning 
the prompt accidental events, the cosmic-induced gammas were observed and the total 
prompt accidental rate is 3.9 /spill/detector, which is 30 times of the proposal 

value. On the other hand, the delayed accidental events rate, of the beam-associated 
gammas, is reduced to 1.2x10“^ /spill/detector, which is 1/39 times of the proposal 
value. This change is due to the reduction of the gammas from the floor by using a 
tighter delayed energy cut and putting 12.5cm thick lead shield below the detector. 
The rate of the beam-associated neutrons decreased to 4x10“^ /spill/detector, which 
is 1/6 times of the proposal value by the on-site measurement. 



Background 

Proposal (/spill) 

New (/spill) 

Prompt 

signals 

e- (i^G(z/e,e-)i^A,.J 

8 . 0 x 10-6 

same 

e- (^^C'(z/a,e-)'^A*) 

3.8x10-6 

same 

e^ (Fg from p~) 

<T(r6 

same 

Cosmic 7 

- 

3.8x10-^ 

Delayed 

signals 

Gamma (Beam associated) 

4.7x10-^ 

1 . 2 x 10-6 

Neutron (Beam associated) 

2.4x10-^ 

4x10-4 

e+ (^^G(z/e,e-)i^A,.,.) 

1 . 2 x 10-6 

1 . 0 x 10-6 

Spallation products 

^T(r^ 

same 


Table 6 : Accidental background summary (per detector). 


The total accidental background is estimated as follows: 


RaccU spill) = X X X 


(4) 


where, i shows some kind of events contributing the accidental backgrounds, 
and are accidental pro mpt a nd delayed rates for each component in Table 

^lifetime lifetime cut (see 3.1). Actually, the AVTX cut efficiency depends on 


each combination of kinds of prompt and delayed like signals, but the AVTX cut 
efficiency for combinations of the z/g prompt signals and delayed like signals of the 
neutron captured on Gd distributes uniformly in the target volume, and it was used 
approximately for all components, which is the same condition as the proposal(2.3%). 
Finally, total accidental background rate (Race) is 6.3x10“'^ (/hour/detector). 


3.2.3 Grand Summary 

Table shows the grand summary of the number of background events at P56. 
The table includes the proposal case with a 17 m baseline and 4000hx4years mea¬ 
surement; the case with new cuts and same operation period but 24 m of the baseline 
(Newl); and the 1.56 times longer period case (5000hx5years), which is the cur¬ 
rent default condition (New2). Values in the Newl case for Fe background from 
/i“(Nj 7 j, ^^G(z/e, background(Ne+) and cosmic-induced fast neutron(N„) 
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are estimated as follows: 


N-p, = 377(proposal) x esL x ^ratio, (5) 

Ne+ = 3S{proposal) x esL x eEd x eAvrx X ^lifetimey (6) 

Nn = 42{proposal) x x (7) 


where, e^L is the squared baseline ratio(17^/24^), eratio is total cut efficiency ratio for 
the signal between cases of this analysis and the proposal (0.81, see |3.2.1 ). eEd and 
^ufetime efficiencies of the tightening delayed energy cut (0.85) and the lifetime 
cut (0.65) for the ^^C(z/e, e“)^^iVg.s, background, respectively (assuming the same case 
as the accidental). e^vTX is AVTX cut between on-bunch and delayed events(0.98). 
^ufetime lifetime cut for the cosmic-induced neutrons (0.65), is total cut 

efficiency ratio for the signal except for the life time cut between cases of this analysis 
and the proposal (0.88, see 3.2.1). 

A clear excess of the Michel electron events from beam-associated neutrons com¬ 
paring the statistical uncertainty was not observed in the on-site measurement. So the 
only upper limit was estimated from the statistical uncertainty (l.lx 10“®/spill/9/rs/detector 
after the on-bunch energy cut). The upper limit is 52.5 times of that in the proposal 
value. Considering the remaining cut efficiencies, except for the on-bunch energy cut 
and the Atprompt cut (total~0.54), the upper limit with a total of 50 tons (2 detectors) 
for 4000 hoursx4 years measurement is 8.2 comparing with 0.2 after considering the 
remaining cut efficiencies (total~0.7) in proposal (the remaining cut efficiencies were 
not considered for the proposal value of 0.3). 

Finally, the dominant background is Fg from (same as the proposal case), 
which is 237, while the signal is 342 in case of the best £t values of LSND, with 
50 tons detector and current default measurement period(5000 hx5 years). 



Contents 

Proposal 

(17m) 

4000hx4y 

Newl 

(24m) 

4000hx4y 

New2 

(24m) 

5000hx5y 

Signal 

sin‘^2e = 3.0 X 10-^ 
(Am^) 

(Best £t values of MLF) 

811 

(3.061/2) 

307 

(2.5el/2) 

480 

(2.5el/2) 

sin^2e = 3.0 X 10"^ 
Am2 = 1.2eV^ 

(Best £t values of LSND) 

337 

219 

342 

background 

Fg from /u~ 

377 

152 

237 


38 

10 

16 

beam-associated fast n 

0.2(0.3x0.7) 

<8 

<13 

Cosmic-induced fast n 

42 

24 

37 

Total accidental events 

37 

20 

32 


Table 7: Grand summary (50 tons in total - 2 detectors). 
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Sensitivity using the Latest P56 configuration 


4.1 Summary of Changed Points from the P56 proposal 

The following points are changed from the P56 proposal: 

• The amount of the accidental background and of the background from the beam 
fast neutrons become realistic due to the MLF 2013BU1301 test experiment; 

• The baseline of the detector is longer than that in the proposal since the detctor 
location should have good background rate to perform the mearurement. In this 
report, we assume 24 meters for the baseline; 

• There are some additional lead blocks under the 25 tons detectors and selection 
cuts to reduce the accidental background while keeping the signal as good as 
possible; 

• Experimental period is expanded to 5 years from 4 years written in the proposal. 
Also the default operation time at the MLF was found to be 5000 hours / year 
during the test experiment, therefore we changed it to 5000 hours / year from 
4000 hours / year written in the proposal. 

For the sensitivity study, the number of events shown in Table are used unless 
it is noted. 


4.2 Methodology of the Fit and Sensitivity 


The methodology of the ht is exactly the same as that of the P56 proposal. It 
uses the profiling method, which treats systematic uncertainties for the normalization 
factors of the signal and the backgrounds correctly. 

In this subsection, we neglect the contributions from the accidental background 
and the background from the beam fast neutrons. The amount of the accidental 
background is about 13.2% of the dominant background (z/g from /r”) and this is 
smaller than that in the proposal. The 90% C.L. upper limit of the background from 


beam fast neutron is less than 5% of the dominant background. Figure 26 shows the 
energy spectra after the selection to be used for the fitting. 

Figure shows the 3a and 5 ct sensitivities with this condition. We find that; 

(1) We will conhrm or refute the LSND region with 5 a (blue) for Am^ > 2eV‘^ region, 
(same as proposal). 

(2) The detector baseline is longer than the proposal, thus the coverage region of the 
3a (green) is larger than that in the proposal in the lower Am^ region. Most of LSND 
90% C.L region can be discussed with 3a. 


Figure 28 shows the variation of the 5cr sensitivity at high Am^(Am^ = lOOeF^) 
region as a function of the operation year (MW x years). This is a good indicator 
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Figure 26: The energy spectra with the signal and the dominant backgrund. Left; 
oscillatation with (Am^, sin^ 26*) = (2.5, 0.003) (MLF best Am^ sensitivity case), 
right; oscillation with (1.2, 0.003) (LSND best fit case). 


whether the P56 experiment can cover whole LSND region above 2 eV^ region, (at 
Am^ = lOOeD^, the left edge of the LSND allowed region by 90%C.L is 0.0017). 

With four years of operation, the sensitivity of the P56 starts to cover the whole 
LSND region with Am^ > 2.0 eV^. 


5 Milestone for the Experiment 

5.1 Accidental Background 


5.1.1 Discussion 


In this subsection, we consider the sensitivity variation with changing the amount 
of the accidental background. This consideration reminds us the importance of the 
background and how to manage them although the current best estimation of the 
background is low enough for the experiment. 

As explained in the section 2.4.1, the accidental background provides an expo¬ 
nential energy spectrum for the prompt energy. If the total amount of the accidental 
background were same as the dominant background (z/g from fi~), the energy spec¬ 
trum of the prompt activities are changed to those in Fig. 
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This has a large impact for the sensitivity since the separation power between the 
oscillated signals and the dominant background is coming mainly from the difference 
of the energy spectra around 50 MeV. The statistical fluctuation of the accidental 
background reduces the separation power even though the energy shape and the rate 


is well determined by beam-off data (less than 1%). Fig. 30 shows the sensitivity at 
the case that there were same accidental background rate as the dominant background 
in 20 < < 60 MeV region. The 5 ct and 3a sensitivities are worsened by about 50%. 
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Figure 27: Sensitivity of the P56 experiment with the latest conhguration (MWx5 
years). Blue line shows the 5a, while the green one shows 3a. The exclusion line of 
the OPERA experiment are also put [5]. The righthand side region of the line are 
excluded with 90% C.L.. 

5.1.2 Conclusion 

The energy shape and rate of the accidental background is found with the 2013BU1013 
test experiment for the hrst time. Thus, taking real data (not only MC) is cru¬ 
cial to know the properties of the background deeper. Also the reduction of the 
background using the hardware and the software including the analysis im¬ 
provement is essential for the experiment if there are larger amount of the accidental 
background than expected. 

During the detector construction phase, we build the two detectors sequentially 
(not in parallel). Thus, we can check the background rate and the energy shape 
with one detector at hrst. During this phase, the best hardware conhguration, e.g. 
shielding, can be also examined. This provides good feed-backs to the other detector 
conhguration and the whole experiment within the short time scale. 


5.2 Experimental Feasibility 

The design of the detector tank is already discussed in the last part of this docu¬ 
mentation. It provides not only the static strength of the tank, but also the endurance 
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Stat.[year] 

Figure 28: The variation of the 5a sensitivity of the P56 experiment as a function 
of operation year (MW x years). 

of the tank against earthquakes. 

We will also prepare a document for other hardware feasibilities. PMTs, electron¬ 
ics, maintenance of the detector will be discussed in the document in detail although 
part of them were already discussed in the proposal. The detector design is similar 
to the reactor experiments (well established detector), therefore the crucial point to 
discuss is the location of the detector and the maintenance of the detector. Note that 
we have to bring the detector out from the MLF building once per year for the MLF 
maintenance at least. 

One detector scheme during the construction phase also helps the realization of 
the maintenance of the detector. 


6 Requests to PAC 

We request the following to PAC; 
(1) Stage-1 approval 
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Figure 29: Energy spectra with the accidental background with the typical oscilla¬ 
tion signals. Left; oscillatation with (Am^, sin^ 29) = (2.5, 0.003) (MLF best Am^ 
sensitivity case), right; oscillation with (1.2, 0.003) (LSND best ht case). Crosses 
correspond to the total, the blue histograms correspond to the dominant background 
(z7e from fi~), and red histograms correspond to the accidental backgrounds assuming 
the rate are same as those of the dominant background in 20 < E < 60MeV region. 



Figure 30: Sensitivity of the P56 with the accidental background (in the the case 
that there were same amount as those of intrinsic z/g-) 
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8 Executive Summary of Background 

In order to understand the backgound rates for the real experimental condition 
at a glance, a summary of the selection critera and their efficiencies at the 25 tons 
detector is shown in this section. Table|^shows the executive summary. Each number 
or selection cut are discribed in the noted sections, and if there are no notes in the 
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table, they are described in the proposal. (Ref [T]) 


Selection 

Signal 

Beam Fast n 

Accidental 

(7+7) 

Accidental 

(7+11) 

1< Atp < 10/is 

74% 

<1.1x10-5 

/9/is/MW/25t 

(90%CL) 

(sec. 2.3) 

3.9x10-^ 
/9/is/MW/25t 
(sec. 2.4.1) 

3.9x10-^ 
/9ps/MW/25t 
(sec. 2.4.1) 

20< Ep <60MeV 

92% 

Atp < Atd < 100/is 

93% 

93% 

1.2x10-5 

/spill/MW/25t 
(sec. 2.4.2) 

1.6x10-^ 

/spill/MW/25t 
(sec. 2.4.2) 

7<Ed< 

(sec. 

C12M 

3.1.1 

[eV 

) 

71% 

71% 

AVTXp rf <60cm 

96% 

96% 

2.3% 

2.3% 

AVTXob 
&& Ey^Xi 
(sec. 

-d <] 

3 <0. 

2.4.2 

10cm 

5MeV 

) 

98% 

98% 

98% 

2.6% 

lifetime 

(sec. 

(L<11) 

91% 

91% 

44% 

65% 

Additional 
(described in 

Ref. PP) 


1/2000 

OB-d spatial(O.Ol) 
#of delayed (0.05) 



Total 

(1 detector) 

38% 

3.1x10“^ 

/spill/MW/25t 

4.6x10-*^ 

/spill/MW/25t 

2.4x10-^ 

/spill/MW/25t 

Events/ 
(2 det 
(sec. 

5year 

ector 

3.2.3 

s/50t 

3) 

) 


< 13 

21 

11 


Table 8: An executive summary of the selection criteria, their efficiencies and pre¬ 
dicted background rates of the 25 tons detector(s ). (Th e final row shows the number 
of the events in 2 detectors as shown in Section. 3.2.3[ Except for the line, all rates 
corresponds to one detector) For the accidental backgrounds, events in the prompt 
region are mainly coming from gammas induced by cosmic rays, while the beam cre¬ 
ates events for delayed region. Note; “OB” means “on-bunch”, “p” corresponds to 
“prompt” and “d” corresponds to the “delayed”. Unless noted, rates and efficiencies 
are discussed in the proposal [T]. 


Note that the accidental background rate is estimated with the following equation; 

RaccUspill) = X X X (8) 

where R is the probability (^backgrounds/spill), and and are the effi¬ 

ciencies of the the AVTXp_rf cut and new additional cuts (the AVTXos-d and the 
lifetime correlation cuts), respectively. Figureshows the distance between vertex 
of the prompt hit and that of delayed hit (AVTXp_rf), which is the reprint from the 
proposal. Using this simple equation, the accidental background rate is calculated 
and is smaller than the intrinsic z7g background. 
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Figure 31: AVTXp_d distributions of i7^ z7g oscillation events (red line) and the 
accidental events (blue line). The gray doted line shows the AVTXp.^ cut condition 
(60 cm). 

A Setup and Calibration of the 500 kg detector 


A.l Setup 

Figure |3^ shows a schematic view of the 500 kg plastic scintillator counters placed 
at the third floor of MLF. It consists of two types of scintillators: 12 pieces of 11.7/13.7 
(trapezoid) x7.6 x 182 cm^ scintillator (ID) and 12 pieces of 16.9/18.8 (trapezoid) 
x7.6 X 182 cm^ scintillator (3D). Each end of the scintillators was viewed by two 
PMTs. Signals from each PMT were recorded by FINESSE 500 MHz FADC. The 
500 kg detector was surrounded by two layers of charged vetoing system, the Inner 
and Outer vetoes. The Inner veto covers the surfaces of the top, bottom and both 
sides of the 500 kg detector. The thickness of the plastic scintillator for Inner veto is 
4.5 cm. The Outer veto surrounds the 500 kg detector and Inner veto, with mostly 2 
layers of 6-8mm thick plastic scintillators. PMT signals from the veto counters were 
recorded by FINESSE 65 MHz FADC with 50 ns RC-£lter. 


A. 2 Calibration 

We used cosmic muons to measure the attenuation length of the scintillator, to 
calibrate the energy and timing. We prepared 4 pairs of scintillator counters to 
trigger the cosmic rays. Figure [33] shows a schematic view of the cosmic muon trigger 
counters. 

To measure the attenuation length of the scintillators, we made some dedicated 
runs in which we changed the position of the trigger counters. Figure M show the 
measured typical attenuation curves for each scintillator type. We measured the curve 
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Figure 32: Schematic view of the 500 kg plastic scintillator detector at the third floor 
of MLF(left: front view, right: side view). The 500 kg plastic scintillators (yellow) 
were surrounded by two layers of charged vetoing system, Inner veto(red) and Outer 
veto (green). 



Figure 33: Schematic view of the cosmic muon trigger counters. We prepared 4 pairs 
of plastic scintillators to trigger cosmic muons. 

























































































and parameterized one for each scintillator. By considering the attenuation length, 
the reconstructed charge was independent from the incident position as shown in 
Fig. 



Figure 34: Typical measured attenuation curve for each scintillator type (ID and 
3D). We measured the curves and parameterized one for each scintillator (magenta 
line). 



reconstructed charge(energy) [MIP] 


Figure 35: Reconstructed charge for different incident positions by considering the 
measured attenuation length. The colors correspond to the events triggered by the 


cosmic muon counters shown in Fig. 33 


The gain of each PMT was also calibrated by using cosmic ray events. Figure [36] 
shows the relative gains of certain PMTs as a function of date. Figure [37] shows the 
relative energy of a certain scintillatoij^with the gain correction as a function of date 
for muon triggered events which are different from the events used for gain correction. 


One scintillator is viewed by four PMTs in total as described in Appendix 


A.l 
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Figure 36: Examples of gain histories for relatively stable(left) and unstable(right) 
PMTs. The black dots show the relative gains for each period, and the red dots show 
the averaged values over stable periods to be used for the gain correction. 


Output from each scintillator is stable within 1-2% with the gain correction. 



Figure 37: Example of the relative energy of a certain scintillator with the gain 
correction as a function of date. 


We also aligned the timing of each PMT. Time offsets were determined to min¬ 
imize the time difference between each pair of PMTs on each end. The velocity of 
cosmic muons passing through the detector was considered. Figure [3^ shows the time 
difference between scintillators for cosmic ray events. The timing of each PMT and 
scintillator was well calibrated. 


A.3 Resolutions 

We evaluated the timing, position and energy resolutions of the 500 kg detector. 
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Figure 38: Time differences between scintillators of cosmic ray events for different 
incident positions. The colors correspond to the events triggered by the cosmic muon 


counters shown in Fig. 33 


We hrst evaluated the timing resolution at each side of the scintillators as a function 
of the equivalent charge. The timing resolution was evaluated by using the timing 
difference between two vertically neighboring scintillators. The total charge difference 
between the two scintillators was required to be less than 20 % of itself, and we can 
thus simply divide the width of the timing difference between those two scintillators 
by \/2 to get the timing resolution. Figure 
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show the typical timing resolutions as 
a function of equivalent charge for two types of scintillators. The obtained timing 
resolution was parameterized with the function: 


®Pi, 


Vc 


(9) 


where C is the normalized charge with the MIP energy at the far end from the side. 
We obtained po = 0.173 ns, pi = 0.246 ns for ID, and po = 0.164 ns, pi = 0.418 ns for 
3D scintillator (typical values). The obtained timing resolutions are mainly limited 
by jitters between channel^ The timing of the scintillator was dehned by averaging 
the hit time at both sides of the scintillator. The typical timing resolutions of the 
scintillators are 0.19 ns for ID and 0.26 ns for 3D at the middle of the scintillators 
for MIP energy. 

We also evaluated the obtained position resolution. The hit position in a scin¬ 
tillator along z-direction was calculated by using the time difference between PMTs 
on both ends of the scintillator. The light velocity in each scintillator was measured 
by using different cosmic trigger events. Figure shows the time difference as a 
function of the hit position. The typical velocity was 14.3 cm/ns. Figure 41 shows 


^ Because we used its internal clock for each FADC, each FINESSE card worked asynchronously. 
We added a reference signal at each waveforms end for each channel every event to align and 
synchronize the FADC cards. 
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charge [MIP@far] 


Figure 39: Obtained timing resolutions as a function of equivalent charge for two types 
of scintillators. The obtained timing resolution was parameterized with Function]^ 


the reconstructed 2;-position for different cosmic trigger events, and Fig. 42 shows the 
difference of the reconstructed z-position between vertically neighboring scintillators 
for MIP energy. By dividing the width by -\/2, we obtained the position resolution 
= 2.6 cm for MIP energy. 


a 



Figure 40: Correlation between 2;-position and the time difference between both ends 
of PMTs. The red line is the fitted function and the slope corresponds to the light 
velocity in the scintillator, 14.3 cm/ns. 


We evaluated the charge (energy) resolution by smearing the output of the pure 
Monte Carlo simulation, and comparing it with data. Figure [43] shows a schematic 
view of the estimation procedure of the charge resolution. We parameterized the 
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Figure 41: Reconstructed 2;-position for different cosmic triggers. The colors corre¬ 


spond to the events triggered by the cosmic muon counters shown in Fig. 33 
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Figure 42: Difference of the reconstructed ^-position between vertically neighboring 
scintillators for MIP energy. By dividing the width by \/2, we obtained the position 
resolution cr^ = 2.6 cm for MIP energy. 


43 



































































charge resolution as follows: 


C 


Po 


Ih 

C 


©P2, 


( 10 ) 


where C is the normalized charge with the MIP energy at the far end from the 
side, po represents the photo-statistics term, pi represents the noise contribution 
and p 2 represents the calibration precision. The charge distribution at 5 different 
positions along the 2;-axis were compared for each PMT. Because only one PMT was 
in interest at a time, we set P 2 = 0- We fitted the charge distributions by changing 
the remaining of parameters at the same time. Light attenuation in the scintillator 
was also considered. The typical energy resolutions of the scintillators are 3.3% for 
ID and 4.5% for 3D at the middle of the scintillators for MIP energy. 



Figure 43: Schematic view of the estimation procedure of the charge (energy) resolu¬ 
tion. The charge resolution was parameterized with Function and five distribu¬ 
tions were htted with the function at once. 


A.4 Detector simulation 


The resolutions described above and other detector responses, such as Birks’ quench¬ 
ing, light attenuation and the threshold effect, were implemented to Geant4 based 
Monte Carlo simulation. Figure 44 shows the deposited energy distribution for a 
scintillator during a beam period with the Monte Carlo estimation overlaid. The 
measured distribution is well reproduced by the Monte Carlo simulation. 
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Figure 44: Energy distribution for a scintillator during beam period with the Monte 
Carlo estimation overlaid. 


A.5 Veto Efficiency 

In order to check the effect of having the Inner Veto (IV) and Outer Veto (OV) 
systems on the MLF third floor background measurement, their particle tagging effi¬ 
ciency was measured. The veto efficiency (e:) is defined as follows: 


Ely 

^OV 


Nb of coincident event OV, IV and Target (Triple Coincidence) 
Nb of coincident event OV and Target (Double Coincidence) 


( 11 ) 


Nb of coincident event OV, IV and Target (Triple Coincidence) 
Nb of coincident event IV and Target (Double Coincidence) 


( 12 ) 


The time window selection for the events was 2.6/rs < Hit Time < 5.3/rs, and the 
energy spectra for different veto condition are displayed in hgure |45j Finally, the veto 
efficiency for different energy ranges is summarized in table where it can be seen 
that a total veto efficiency is bigger than 99.5%. In these table and figure, it is also 
possible to see that muons will deposit an amount of energy larger than 60 MeV and 
almost of the totality of them will thus be vetoed. 


B Small size detectors to measure the PID, Rate 
and Energy of Prompt BKG 


In a neutrino oscillation experiment, such as one looking for flavour appearance, 
it is important to understand the backgrounds that can falsify a true neutrino in¬ 
teraction. Even though the selection cuts are applied to maximize the Inverse Beta 
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Figure 45: Target energy spectra for different veto conditions. 


Table 9: Summary of IV and OV efficiency for different energy ranges. 


Energy Range [MeV] 

£iv 

^ov 

(IV or OV) 

20 < E < 60 

96.8±0.2% 

94.1±0.2% 

-99.8% 

60 < E < 100 

99.5±0.04% 

96.2±0.1% 

—99.9% 

100 < E <140 

99.6±0.07% 

95.1±0.3% 

—99.9% 
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Decay (IBD) sample in the data, there are accidental coincidences that can pass these 
cuts. Moreover, an energetic neutron entering the detector, hereafter called fast-n, 
can create a correlated prompt and delayed event pair when it scatters a proton and 
it is thermalized by the liquid scintillator, hnally being captured by a Gadolinium 
nucleus, mimicking an IBD event. Natural occurring high energy gammas and neu¬ 
trons are secondary products of the cosmic rays interactions in the upper atmosphere, 
from their hadronic and electromagnetic components. In addition, a fast-n can in¬ 
teract with atomic nucleus of the materials in the detector surroundings, leaving the 
nucleus in an excited state that will produce a high energy gamma during its de¬ 
excitation; or the nucleus can absorb the neutron, which will also emit gammas in 
the process. Therefore, the sources of accidentals and correlated background must be 
understood and their rates predicted. This section covers the use of small detectors 
to study such backgrounds. 


B.l Measurements at Tohoku University 

In the laboratories of the Research Center for Neutrino Science of Tohoku University 
(RCNS), several measurements were performed and their data compared with Monte 
Carlo simulations (MC) in order to understand those backgrounds. Later, as it is 
described in the next section, these measurements results were compared with the 
ones performed at the third floor of MLF facility and extrapolated to the detector 
type and size being proposed in this document. 


B.1.1 Nal 

A Sodium Iodide (Nal) scintillator counter (a cylinder with 2” of diameter and 2” 
height) was hrst used to measure the high energy gamma ray flux that contributes 
to the prompt event of the accidental background. The Nal was fully surrounded by 
plastic scintillator counters (six identical rectangular cuboids of 20 x 90 x 5 cm) acting 
as a veto, in order to ensure the purity of this neutral component. Figure [46] shows the 
data used to calibrate the MC for this set-up, while the right plot of hgure [47| shows 
the full energy range and the MC with its components. After the environmental 
gammas, cosmic muons and neutrons were generated, a missing component of the 
spectrum was identihed to be high energy gammas, which energy spectrum can be 
described by two exponential functions with decay constants of 3 and 26 MeV. These 
components rates was measured to be 150 and 25 Hz/m^, respectively. 

Finally, the effect of adding the veto system was studied. The left plot of hgure [47| 
shows the data veto spectrum compared with total the MC and its components. 
The doted line represents the cut value used to reject events depositing considerable 
(bigger than 500 default units of charge [DUQ]) amounts of energy in any part of the 
veto. This plot also shows that there is an excess when comparing the MC and data 
histograms between the cut value and 5000 DUQ. One possible explanation for this 
excess is that it is caused by Michel-e“, and that the simulation does not apply well 
the low energy muons that decay and produce such events. Thus, the 500 low value 
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Figure 46: Energy spectra of the low energy (< 3 MeV) environmental gammas 
using a Nal counter, where the potassium and thallium peaks were used for the MC 
calibration. 


I Veto Charge | 


I Nal Energy | 



Charge 


—•— Data 

Iv a; y-\ Env.Gamma{MC) 

I I Cosmic Muon(MC) 

I I Cosmic Gamma(MC) 



Figure 47: Full energy range spectra for Nal counter and its veto system. The 
left plot shows the charge distribution for data and MC of the veto system where 
the doted line represents the veto cut value, while the right plot shows the full Nal 
energy spectrum where the MC is decomposed in its components: two exponentials 
with decay constants of 3 and 26 MeV for gammas and cosmic muons. 
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Nal Energy 



Figure 48: Nal energy spectrum for the remaining events after applying the veto 
cut. 


cut was used, to avoid this unknown region. The remaining events, after applying the 
veto, are shown in hgure where it can be seen that on the prompt energy window 
(20 to 60 MeV) the events are mainly gammas. 


B.1.2 NE213 


The next important component of the natural background that was studied is the 
fast-n. Thus, a liquid scitillator with Particle Identihcation (PID) capabilities was 
chosen. A cylindrical aluminium housing, with white painted inner walls, of 613.6 
cm^ was hlled with NE213, closed with a glass plate and attached to a 5 inches 
PMT (R1250-03). This detector where surrounded by the same plastic scintillators, 
described in the previous subsection, making the veto system. The detector and its 


veto are shown in hgure 49 


The NE213 is a well known scintillator for Pulse Shape Discrimination (PSD). 
The ratio of its fast and slow time decay constants is dependent on the type of the 
incident particle. Therefore, it was used to separate the neutron (recoiled proton) 
from gamma (mainly Compton electrons and positrons from pair production) signals. 
As for the Nal case, the NE213 set-up was compared with MC for tuning its 


parameters. On hgure 50 the photon-electron (PE) spectra of data and MC for envi¬ 
ronmental gammas and a Cobalt-60 (®°Co) source are presented. While the nominal 
value is 13565 PE/MeV the plots shows that a tuned MC agrees well for a value of 650 
PE/MeV. This diherence is mainly due the Quantum Efficiency (QE) of the PMT, 
the rehection inefficiency of the walls and leakage through the optical contacts. 

Another important parameter tuned was the scintillator Birks constant. A well 
dehned Birks value is necessary due the proton elastic scattering by neutrons, other¬ 


wise its spectrum would have an incorrect shape. Figure M shows the neutron-like 
and gamma-like events for data and MC when using an Americium-Beryllium (AmBe) 
source, where the MC value for the Birks constant is 0.100 mm/MeV, while a refer- 


49 




















































Figure 49: Scheme of the NE213 detector surrounded by the veto system formed by 
plastic scintillator blocks. 
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Figure 50: Energy spectra of the environmental gammas and ®°Co using the NE213 
counter. 
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enced value is 0.107 mm/MeV [B]. The AmBe source is a well used source for neutron 


MC n-like spectrum 



MC 7-like spectrum 



Figure 51: AmBe data and MC comparison. The left plot is the events tagged as 
neutrons and the right plot for the ones tagged as gammas. The agreement between 
the spectra for the neutrons shows that a good value for the Birks constant was 
chosen. 

studies. When an Am nucleus emits an alpha particle, it can interact with the Be 
producing free neutrons and gammas. The method to separate the two components 
of the AmBe spectrum, was to divide the tail integral of the electric pulse by the total 
integral. Since the slow decay component of the scitillator has a higher contribution 
for neutrons than for gammas, this ratio variable {Qtaii/Qtot) has a higher value as 
well. 

As shown in hgure the neutron component is distinguishable from the gamma 
one, when plotting the Qtail/Qtot versus the event reconstructed energy. Although 
the shape of the data and MC scatter plots does not agree well, due to difficulties 
on calibrating the MC for timing and pulse creation, both particles can be selected 
when applying a correct energy and BID cut. 



Figure 52: AmBe scatter plot of Qtaii/Qtot versus the event reconstructed energy. 
The left plot shows the data and the right one the MC. The doted horizontal lines 
represent the cut values to separate the neutron and gamma components. 


The next step was to calculate the efficiency of the data sample’s PSD cut. The 
data sample in the energy range (from 7 to 61 MeV) was divided in six regions: 1) 7 
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to 11, 2) 11 to 16, 3) 16 to 23, 4) 23 to 30, 5) 30 to 38, and 6) 38 to 61 MeV. For each 
region an one dimension histogram was constructed for the PSD variable and htted 
by a two peak Gaussian function. Each peak of the Gaussian function represents 
the neutron and gamma components. Then, a cut value was chosen to distinguish 
both components and its efficiency was calculated based on the htted Gaussian. This 
method is depicted by hgure where the histogram for each region is displayed 
together with the htted Gaussian and the neutron selection cut value represented by 
a vertical line. The neutron and gamma component of the Gaussian is shown in blue 
and red respectively. 


Region 1 Region 2 




Region 3 



Region 4 Region 5 Region 6 





Figure 53: Qtaii/Qtot histograms for each region dehned in the text. The data 
was htted by a two peak Gaussian. The red and blue lines represent the gamma 
and neutron component respectively. The vertical line is the cut value to select the 
neutron events from the data for each region. 


In hgure 54 the Qtail/Qtot cut value, cut efficiency for neutrons and gamma con¬ 
tamination at the neutron sample, for each region are shown. 

Finally, as performed with the Nal, hgure 55 shows the veto charge distribution 
on the left and NE213 reconstructed energy on the right, where the data and MG 
are presented. On hgure 56 the vetoed selected neutron and gamma components, 
as explained above, can also be seen. On this hgure it is possible to see that the 
components agree well between the data and the MG model. 

With this set-up, the fast-n hux detected above 20 MeV is of 1.28±0.05 x 10“^ Hz 
(statistic uncertainty), while the MG gives a rate of 1.12 x 10“^ Hz. These numbers 
are close to the one used in the hrst proposal. For gammas, in the same energy range, 
the measured rate is of 1.18 ± 0.04 x 10“^ Hz (statistic uncertainty), while the MG 
gives 0.95 x 10“^ Hz. Therefore, the data and MG above 20 MeV for both gammas 
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Figure 54: Qtail/Qtot cut value (left), cut efficiency for neutrons selection (center) 
and gamma contamination at the neutron sample (right), for each regions defined in 
the text. The cut efficiency for neutron events selection is dehned as the ratio between 
the events above the cut value and the total integral of the fitted Gaussian function. 
Similarly, the gamma contamination of the neutron sample is the ration between the 
gamma events above the neutron cut and the selected neutrons. 
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Figure 55: Data and MC comparison for the NE213 set-up. The left plot shows the 
veto charge distribution and the right plots shows the NE213 reconstructed energy 
distribution. 
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Figure 56: Reconstructed energy distribution for events with low or none energy 
deposition in the veto. The data and MC components of neutrons and gammas, 
selected as explained in the text, are also compared. Data events without veto applied 
are also shown for comparison. 


and fast-n agrees within 20% of the uncertainty. For the MC presented here and 
for the one in the hrst proposal, the same cosmic fast-n generator was used. The 
generator’s flux and spectrum were tuned with Tohoku University’s Reactor Monitor 
detector, composed of 200 litters of liquid scitillator. 

C Gamma Ray Measurements with a Small Plas¬ 
tic Scintillator 


C.l Setup 

To validate the MC model for gamma background, we measured the event rate at 
Point 2 with a small plastic scintillator. The scintillator dimenssions are 90 x 21 x 4.5 
cm^ and a PMT is attached to side of the scintillator with grease as shown in Fig. 
The PMT waveforms were digitized with a 6 /rs time window and recorded through 
a FINESSE 65 MHz FADC. The trigger system was same as the 500 kg detector. 
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C.2 Lead Shield Configurations 

We measured the number of events induced by gammas with 4 shield conhgura- 
tions as shown in Fig. 


1. Without lead shield (top left of Fig. 58, labeled “Without lead”) 

2. 5-cm thick lead shield surrounding nearly all over the scintillator (top right of 
Fig. 58, labeled “Lead dvr”) 
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PMT attached 
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Figure 57: Schematic of the small plastic scintillator 


3. 5-cm thick lead shield under the scintillator (bottom left of Fig. 58, labeled 
“5-cm thick lead”) 


4. 10-cm thick lead shield under the scintillator (bottom right of Fig. 58, labeled 
“10-cm thick lead”) 


C.3 Energy Calibration 

For the energy calibration, we used the peak of cosmic muons. Using MC, we es¬ 
timated that the peak energy is about 9 MeV. We that the assumed energy resolution 
around 9 MeV for MC is 18 % mainly caused by attenuation of light. A calibrated 
energy spectrum taken while the accelerator was off and a MC energy spectrum made 
by cosmic ray muons with energy resolution are shown in Fig. [59 


C.4 Measurements 


C.4.1 Energy spectra of 4 configurations 


As shown in Fig. 60, data includes beam bunch timing. To avoid effects of beam 
activity, we use data 3.5 /is after the gate start time. Converting FADC count into 
energy, the energy spectra of 4 conhgurations were made. To remove the effect of 
environmental gammas and cosmic rays, we subtracted the accelerator off data from 


these spectra. The energy spectra after the subtraction are shown in Fig. 61 
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Without lead Lead~47r 

live time: 12.5hour live time: 24.7hour 



5cm-thick lead 
live time: 2.4hour 



1 Ocm-thick lead 
live time: 40.8hour 



Figure 58: Shield configurations of the small scintillator. 



Figure 59: Energy calibration using cosmic ray muon. The black line shows the 
measurement and the red line shows MC. 
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Figure 60: Example of FADC count vs time plot (’’Without lead” configuration). 
Around 1-2/is, beam bunch structure is clearly shown. To avoid effects of beam ac¬ 
tivity, we use data 3.5 /rs after the gate start timing. FADC count of 150 is equivalent 
to ~10 MeV. 
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Figure 61: Energy spectra of 4 conhgurations after subtraction of accelerator off data. 
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Figure 62: Comparison between measurements and MC. The legend numbers are 
the rejection powers. 


C.4.2 Validation of the MC model: rejection power and energy spectra 

For each configuration, the energy spectrum was reproduced by the MC. Accord¬ 
ing to the MC model, gammas are generated uniformly in a 5 x 5m^ fioor surface, 
where the small scintillator was centred. Deposit energy was smeared with a 18% 
energy resolution. 

To validate the MC model for gamma background, we compared rejection powers 
of measurements and MC. The definition of rejection power, r, is; 


r = 


N^ 


Without lead 




config 


(13) 


where Awithoutlead is the event rate of 4 < F'[MeV] < 10 at “Without lead” configu¬ 
ration and Aconfig is that at other configurations (“Lead~ dvr”, “5cm-thick lead”, 
“lOcm-thick lead”). In addition, we checked the energy spectra. The rejection 
powers and energy spectra comparing between measurements and MC are shown 
in Fig. 1^ For normalization, the energy spectra of MC were multiplied by a fac¬ 
tor, A“utiead/^Widhoutieadi ^here A“utiead is evcut rate of measurement and 
^Without lead i^ event number of MC in energy region of 4 < A [MeV] < 10 at “Without 
lead” configuration. 

The MC model explains the measurements well because the rejection powers for 3 
configurations are consistent between measurements and MC, and the energy spectra 
shape of MC also reproduce the measurements below 9 MeV. Therefore the MC model 
was validated and the lead shield effect was confirmed. 
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